Global Persistence Exponent in Critical Dynamics: Finite Size induced Crossover. 
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We extend the definition of a global order parameter to the case of a critical system confined 
between two infinite parallel plates separated by a finite distance L. For a quench to the critical 
point we study the persistence property of the global order parameter and show that there is a 
crossover behaviour characterized by a non universal exponent which depends on the ratio of the 
system size to a dynamic length scale. 
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Global persistence exponent for non-equilibrium crit- 
ical dynamics was introduced a decade ago[T|, following 
the emergence of similar exponents in the evolution of 
Ising spins in one and higher dimensions and the 

evolution of a diffusing field ^ from random initial con- 
ditions in different dimensions. The simplest system ex- 
hibiting persistence is the random walk in one dimension 
Q. Since Brownian motion under restrictive geometry 
has been of experimental interest lately 0], the persis- 
tence problem was addressed under those situations Q. 
It was seen that the power law decay for the infinite sys- 
tem acquired an exponential correction for the confined 
system (confinement by walls or harmonic forces). This 
was in contrast to the finite persistence probability ob- 
served by Manoj and Ray for finite size systems exhibit- 
ing critical dynamics. The quench carried out by Manoj 
and Ray [9] was, however, deep into the ordered region. 
For a D-dimensional Ising model, starting from a random 
initial condition, they quenched the system to T = and 
allowed the spins to evolve according to Glauber dynam- 
ics. Domains began forming and when it happened that 
the domain size became larger than the system size , 
then the persistence probability attained a finite value. 
The global persistence exponent of Majumdar et al [ij 
was defined differently. It referred to the quench from 
a high temperature to T = Tc, the critical point of the 
system and considered the global order parameter. The 
individual spins flip rapidly and the probability of not 
flipping in an interval has an exponential tail. It is only 
when the global order parameter is considered, that one 
flnds the power law tail. In this situation if we consider 
a flnite size system, then for a sufficiently small system 
size (smaller than the appropriate "dynamical" length 
scale), the global order parameter will no longer find it 
so difficult to "overturn" and an exponential tail could 
be expected just as happened with the Brownian motion 
in restrictive geometry. In this note we use the spherical 
limit to establish our result. 

We consider the usual Landau Ginzburg free energy F for 
the A'^-component order parameter (pi {i — 1, 2, ....N}, in 
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a 3-dimcnsional space, that is 
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F = / d^f[-0,0, + -(V,0,)(V,0O + -M^<i}i){(pAl)]■ 
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The corresponding Langevin equation is given by 
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where ^ is a Gaussian white having correlation 
< S,{r,t)i{r',t') >= 2rS{f- ?)6{t - t'). 
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Since we will be using spherical limit, it makes sense to 
work in D = 3 directly. The range of validity of the 
spherical approximation if for 2 < D < 4 and hence 
Z? = 3 is the natural choice. 

The confinement is taken to be in the z-dircction and 
the orthogonal space has two dimensions. The confining 
is in the form of two "parallel plates" at z = and at 
z = L, where Dirichlet boundary conditions hold pjO|- 
. The other two dimensions are infinitely extended. 
The decomposition of (pi (r , t) is now in terms of Fourier 
transform in two dimensions and a Fourier series in the 
z-direction, so that 
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and the linearized Langevin equation becomes 
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i,„{k) -Trpi^nik), (5) 



in the non-interacting limit, w 

2 -^ 

n=l,r=—j^,k — gives us 



0. For the choice of 
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At the critical point for the confined system (r = — tt^/L^ 
represents the mean field expression of the critical point), 
the lowest mode (fc = 0,n = 1) undergoes a Brownian 
motion, corresponding to a persistence exponent 6 ~ 0.5. 
For the finite size system, we identify 0i,i(O) as the global 
order parameter. 
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To work in the spherical hmit we write Eq.([2]) as 



' Ak)-T{r + -N <4>^ >) 



we arrive at 
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= l/[s + -j^ + 4r\{J{0, L) - J{s, L)}], (15) 



(f>^^n{k) + uO{l/N)(P,^n{k) + U^, n, t).(7) where 
Since {N < 0^ > -J2j'l^]) is of 0(1) and hence in the 



hmit N ^ oo (spherical limit), we have, for any i. 
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Defining a{t) — — r(r + u < cj)^ >), we can write 
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The solution for (j}n{k^t) can now be written as 



(j)n{k,t) = e 



ak,t') 



(f)n{k,0)e' 







-r(fc^ 



'-)t+b(t) 



(10) 



where b{t) — dt'a{t'). The long time dynamics is dom- 
inated by the noise containing term and < 0^ > in that 
limit is given by. 
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where 
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The dynamics of g{t) is given by 

g = -2gh{t) = -2ga{t), 
= 2gT{r + u< (fp- >), 

= 2rTg + AuT [ dt'g(i') V 



The critical point is now defined by the zero of the coef- 
ficient of the k — 0,n — 1 component of 4>n{k) in Eq.® 
and thus 



u < (b" >= - 
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If we consider the Eq. p^ at the critical point, then in 
the terms of the Laplace transform 



~9{s) 



g{t)e-''dt, 



with 
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A J = J(0,L) - J(s,L), 
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For L — > 00 we can write 
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The first correction to A J for L — > 00 is given by the 
second term in Eq. (|17p which becomes 



d2fc 



kdk 



(27r)2 (s-f 2r/c2)(2r/c2) 2t: J {s + 2rk^)(2Tk 

dz 



1 

8^ 



In 



^ + s/2r 



(19) 



The divergence at the lower end needs to be cut off and 
this is done by recognizing that the lowest value of k is 
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0{L ^) and to the leading order the integral is 
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Consequently for L finite but much greater than s 
the expression for A J becomes 
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We now need to explore the limit L(^)-'^/^ <C 1. To do 
this we return to Eg. pTp . perform the two dimensional 
k integration and write 
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Since the denominator of Eq. (|15p already contains a term 
linear in s this limit A J will not reveal any additional 
feature. For small value of s, we can now write for 
L(^)i/2 » 1 



2lT 



wrL(s/2r)i/2 



1 



ln(L2s/2r) 



2tt 



2T 



1/2 



L(s/2r)i/2 
2r\ (L^s 



(22) 



The real time behavior is obtained by inverting the 
Laplace Transform of g{s) and we have 
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At this order the expression for a{t) and b{t) becomes 
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We note that for L large enough so that ^ 1, a{t) 
jT, with the first correction given by 
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For -^l <C 1, we can write a(t) as where e(t) is the 
quantity in brackets in Eq. (|25p and is slowly varying func- 
tion in the range considered. 

The global mode (j)i{0) now satisfies the equation of 
motion (see Eq.®, 
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Under the transformation (/«i(0,t) = e""'''/-^'t'(*)/'*V(0 
and making the slowly time varying approximation 
whereby (e/e)i Int is considered significantly smaller than 
unity (that is Tt/L^ reasonably smaller than unity) we 
arrive at 

^{t) = e"^'/^'t-^(*)/4^(t) 
With the transformation of variable t = we get 

^ir)^=e^'^'^'^'t-<'^/'m-f{r) (27) 



The correlation function < /(t)/(t') > will be delta cor- 
related in T-space provided 
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and Eq. (|27|) becomes 
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Since the size dependent correction in e is 0{L^^), we 
can drop the last term to the leading order and write as 
the first effect of the finite size, the relation 
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The persistence probability for the process of Eg. ([291 
goes as T~^/2 and hence in the actual time variable 
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The decay is clearly hastened at a finite value of L. 

What happens is /Tt becomes smaller than unity? 
Returning to Eq. (fT5|) and Eq. ([2T|) . it is now clear that 
the leading behavior of g{t) is e^^'" leading to h{t) — 



^-^t and a{t) = This implies a dynamics 
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The associated p{t) is known from ref. [3] to be 
Pit) = 



L2 



(32) 
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A combination of the forms of Eq. (l3T|) and Eq. p3)) can 

be achieved by 
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where 
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For 3> Tt, we have the result of Majumdar et. ah 
[H, that is p{t) ^ while for ^ Ft, we regain 

Eq.dSH). 



We note that the critical relaxation rate for a system 
governed by Eq. Q goes as Tk^ , where k is the wavenum- 
ber of the fluctuations. For the finite sized system the 
minimum value of k is 0{L~^) and the lowest frequency 
is given by TL^^ . For any arbitrary time scale 't', the 
relaxation rate allows us to define a dynamic length scale 

= The ratio L? /Tt which has featured so promi- 

nently in our discussion is thus the ratio L'^/l'^- The 
results of Majumdar et al [l| are for the limit L ^ Id- 
Our Eq. (|34p is an attempt to capture the entire range 
from L ^ and L ^ 1^. 
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